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Abstract. We investigate the quarkonia correlators in QCD with no light quarks within a potential model 
with different screened potentials. Our results for the temperature dependence of the charmonium and 
bottomonium correlators are qualitatively consistent with existing and preliminary lattice results. We 
identify however, a much reacher structure in the correlators than the one seen on the lattice. 
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1 Introduction 

Matsui and Satz conjectured that color screening pre- 
vents the binding of heavy quarks above the deconfine- 
ment transition [Q. This idea inspired the intense inves- 
tigation of heavy quark bound states in hot strongly in- 
teracting medium, as this can test deconfinement. Refine- 
ment of this idea led to the identification of a hierarchy 
in the dissolution pattern, meaning, that higher excited 
states disappear earlier [2IEI- 

Recently available numerical analysis of quarkonia cor- 
relators and spectral functions carried out on the lattice 
for quenched QCD 4,5,6 provided unexpected results: 
The ground state charmonia, IS J/ip and ?y c , survive at 
least up to 1.5T C . Not only that these states do not melt at, 
or close to T c , as it was expected, but lattice also found 
very little change in their properties when crossing the 
transition temperature. In particular, the mass of these 
states show almost no thermal shift. Furthermore, lattice 
results indicate [S] that properties of the excited IP states, 
Xc and Xc > are seriously modified above the transition tem- 
perature, and that these states are dissolved already at 
1.1T C . Preliminary results from the lattice for the temper- 
ature dependence of the bottomonia correlators are now 
also available and predict the melting of these states 
at much higher temperatures than the charmonia states. 

The above features are in contrast with existing poten- 
tial model studies, which predicted that the J/tp would 
disappear at around 1.1T C (see for example [3]). There 
have been some recent attempts to understand these re- 
sults: In [8] the J/tp was considered as a strongly coupled 
color-Coulomb bound state, and found to survive until 
2.7T C . In |H for the J/tp a spontaneous dissociation tem- 
perature of 2T C has been found using a potential fitted 
for lattice results. In these studies however, only the ques- 
tion of existence or non-existence of quarkonium binding 
was addressed. No comparison of the quarkonium proper- 
ties was made with the available lattice data, although it 



is expected that these are strongly modified close to the 
corresponding melting temperatures. 

Since the calculation on the lattice of meson spec- 
tral functions at finite temperature is difficult even in 
quenched QCD, we study the quarkonium correlators in 
QCD with one heavy quark (either c or b quark) in terms 
of a potential model. This allows us to make a direct con- 
nection between lattice calculations of quarkonium at fi- 
nite temperature and potential models. 



2 The Model 

We investigate the quarkonium correlators in Euclidean 
time at finite temperature, as these are directly calculable 
in lattice QCD. The correlation function for a particular 
mesonic channel H is defined as 



G H (r,T) = (j H (r)jU0)) 



(1) 



Here j H = qr H q , and r H = 1, 7^,75,7^75 corresponds 
to the scalar, vector, pseudo-scalar and axial vector chan- 
nels. Using its relationship to the retarded correlator, the 
following spectral representation for Gh (t, T) can be de- 
rived 



G(t,T) 



cosh [lu (t — T^f) 
sinh [£] 



(2) 



Here we need to specify the spectral function at finite 
temperature. We do this by following the form proposed 
in |10j for the zero temperature spectral function 



(uj 2 - M?) + m uj 2 e (w - s ) . (3) 



The first term contains the pole contributions from reso- 
nances, while the second term is the perturbative contin- 
uum above some threshold. For QCD with light dynamical 
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Table 1. Parameters at T = 



Q 


a 


m c 


m b 


SOc 


SOi, 


0.471 


0.192 GeV 2 


1.32 GeV 


4.746 GeV 


4.5 GeV 


11 GeV 



quarks it is natural to identify the continuum threshold so 
with open charm or beauty threshold. For the case of one 
heavy quark only, which is considered here, the value of 
so is somewhat arbitrary. We choose this threshold as the 
energy above which no individual resonance are observed 
experimentally. These are shown in Table ^ Om finite 
temperature model spectral function has the form given 
in |J3J, with temperature dependent decay constant Fi(T), 
quarkonium mass M,(T), and threshold sq(T). For the pa- 
rameter mo we take the value 3/(87r 2 ) calculated in lead- 
ing order perturbation theory The threshold above 
which quarks propagate freely is given by twice the ef- 
fective quark mass, i.e. temperature-dependent pole mass 
sq(T) = 2m p (T). We consider only the case of zero spa- 
tial momentum, and therefore do not have to take into 
account contributions from Landau damping to the spec- 
tral function. 

The remaining parameters of the spectral function can 
be calculated using a potential model. The bound state 
mass is given by Mi = 2m c ^ + Ei, where Ei is the binding 
energy, and m Ci b are the constituent charm and bottom 
quark masses given in Table ^ At leading order in the 
coupling and inverse mass the decay constant can be re- 
lated to the wave function at the origin for the vector and 
pseudo-scalar channel, and to its derivative for scalar and 
axial- vector channels We obtain the binding energies 
and the wave functions by solving the Schrodinger equa- 
tion. We solve this with a temperature dependent screened 
Cornell potential, first considered in [2] 



V ^ T ) = --e-^ T > + 



(1 



(4) 



We parameterize the temperature dependence of the screen- 
ing mass as /i = (0.24 + 0.31 ■ (T/T c - l))GeV with T c = 
0.270 GeV critical temperature. Here a is the coupling 
and a is the string tension. Their values as well as quark 
masses were obtained in |13] by fitting the zero tempera- 
ture quarkonium spectrum and given in Tabled 

In our analysis we also used a screened potential with 
a quite different functional form. This was obtained by fit- 
ting the quenched lattice QCD data on the internal energy 
for T > 1.18T C We have not considered the internal 
energy closer to T c as it is not clear if it can be identified 
with the potential. 

The potential (J3J is shown in FigureHJfor different tem- 
peratures above deconfincmcnt. A common feature of all 
screened potentials is that they reach a finite temperature- 
dependent asymptotic value, V^T). This contributes an 
extra thermal energy to the quark anti-quark pair in the 
plasma |15j . This thermal internal energy thus leads to an 
effective quark mass given by m p (T) = m Ct b + V^o(T)/2. 
Since T4o(T) is decreasing with temperature, also so(T) = 
2m p (T) will decrease accordingly, the effect of which man- 
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Fig. 1. Cornell potential at T — and screened Cornell po- 
tential for different T > T c . Here a = 0.192 GeV is the zero 
temperature string tension. 



ifests in the temperature dependence of the correlators, 
as we discuss in the following section. Such a decrease 
in the pole mass was observed in lattice calculations |16|. 
when calculating quark and gluon propagators in Coulomb 
gauge. While we assume that above the threshold quarks 
and antiquarks propagate with the temperature depen- 
dent effective mass defined above, quarks inside a singlet 
bound state will not feel the effect of the medium, and thus 
will have the vacuum mass. Therefore in the Schrodinger 
equation we use the zero temperature masses of c and b 
-quarks. Above deconfincmcnt the quarkonium can also 
dissociate via its interaction with gluons ^Z]- This effect 
leads to a finite thermal width, which we will neglect in 
the following analysis. 



3 Results 

Here we concentrate on results obtained in the scalar and 
pseudoscalar channels. Details of the analysis, together 
with the study of the vector channel, where transport ef- 
fects are also considered, will be presented elsewhere jTSJ ■ 

First, let us discuss the properties of the quarkonia in 
the deconfined medium. In Figs. and El we present the 
temperature dependence of the masses and the amplitudes 
F(T) (i. e. the wave function or its derivative at the origin) 
of the different quarkonia states. One can see that while 
the masses do not change substantially with temperature, 
there is quite a strong drop in the corresponding ampli- 
tudes. While the small shift in the quarkonia masses above 
the deconfinement temperature is consistent with lattice 
data, the decrease in the amplitudes is neither confirmed, 
nor ruled out by existing lattice data. 

In order to make a direct comparison with lattice re- 
sults we normalize the correlation function to the so-called 
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Fig. 2. Temperature dependence of quarkonia masses normal- 
ized to their corresponding zero temperature masses. 
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Fig. 3. Temperature dependence of quarkonia amplitudes nor- 
malized to their corresponding zero temperature amplitudes. 



reconstructed correlators [5] 



G recon (T,T)= / dua(u,T = 0)K{T,u,T). (5) 
Jo 

Studying the ratio G/G recon can indicate the modifica- 
tions of the spectral function above T c . Deviations of this 
ratio from one is an indication of medium effects. 

The results for the scalar x° c ancl pseudoscalar rj c char- 
monia correlators are presented in Figure^ There is a very 
large increase of the scalar correlator (upper panel) which 
is in qualitative agreement with lattice data [Sj • The scalar 
correlator above deconfinement is enhanced compared to 
the zero temperature correlator despite the fact that the 
contribution from the \c state becomes negligible. This 
enhancement is due to the thermal shift of the continuum 
threshold. For the pseudo-scalar correlator (lower panel) 
we also find a moderate increase. This increase of the cor- 
relator is again attributed to the decrease of the threshold 
with temperature. This extra feature, i. e. the significant 
contribution to the correlator from the continuum due to 
threshold reduction, is not detected on the lattice for the 
pseudoscalar channel. 
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Fig. 4. Ratio of charmonia correlators to reconstructed corre- 
lators for scalar (upper panel) and pseudoscalar (lower panel) 
channels. 




Fig. 5. Ratio of the pseudoscalar charmonia correlator to the 
reconstructed correlator that includes IS and 2S states. 



A qualitatively similar behavior was obtained for the 
bottomonia states, xt an d Vbi as illustrated in Figure 
The behavior of the scalar bottomonium channel is very 
similar to that of the scalar charmonia, even though con- 
trary to the Xc state the Xb survives until much higher 
temperatures. This is due to the fact that the shifted con- 
tinuum gives the dominant contribution to the scalar cor- 
relator. In the pseudoscalar channel wc see a drop in the 
correlator at large enough compared to the zero temper- 
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ature correlator which is due to the reduction of the am- 
plitude. This is presumably due to the fact that IS bot- 
tomonia are more deeply bound thus there is a more clear 
separation between the pole and continuum contributions 
to the correlator. 
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Fig. 6. Ratio of bottomonia correlators to reconstructed corre- 
lators for scalar (upper panel) and pseudoscalar (lower panel) 
channels. 



In the analysis presented so far we have considered 
only the lowest meson state in a given channel. This is 
because we wanted to make contact with existing lattice 
data in which the excited states in any given channel were 
not yet observed. To study possible effects due to higher 
states in a given channel we have included also the 2S 
state in the pseudoscalar channel in the spectral function 
10. Since the 2S state is melted at T c , we expect a drop in 
G/G recon compared to the previous case, where only IS 
state is included. Figure El shows the temperature depen- 
dence of the charmonium pseudoscalar correlator obtained 
when both the IS and the 2S states are accounted for. As 
expected, in this case we identify a 25 % reduction in the 
correlator, attributed to the melting of the 2S state. This 
has not been seen on the lattice, since data involving the 
2S state are still unavailable. 

To test the robustness of our results we reconsidered 
the analysis with a different parameterization of the screen- 
ing mass. For temperatures T > 1.18T C we also used the 
potential fitted to the internal energy calculated on the 



lattice. Our results proved not to be sensitive to the de- 
tailed form of the potential, suggesting that they are based 
on very general physical arguments |18|. 

In conclusion: We have analyzed mesonic correlators 
using model spectral functions based on the potential pic- 
ture with screening. We have found qualitative agreement 
with lattice data. On the other hand the temperature de- 
pendence of the correlators show a much reacher structure 
than the one seen on the lattice. 
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